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Previous models of the cerebrovascular smooth muscle cell have not addressed the interaction be- 
tween the electrical, chemical and mechanical components of cell function during the development of 
active tension. These models are primarily electrical, biochemical or mechanical in their orientation, 
and do not permit a full exploration of how the smooth muscle responds to electrical or mechanical 
forcing. To address this issue, we have developed a new model that consists of two major compo- 
nents: electrochemical and chemomechanical subsystem models of the smooth muscle cell. Included 
in the electrochemical model are models of the electrophysiological behavior of the cell membrane, 
fluid compartments, Ca^"*" release and uptake by the sarcoplasmic reticulum, and cytosolic Ca^^ 
buffering, particularly by calmodulin. With this subsystem model, we can study the mechanics of 
the production of intracellular Ca?^ transient in response to membrane voltage clamp pulses. The 
chemomechanical model includes models of: (a) the chemical kinetics of myosin phosphorylation, 
and the formation of phosphorylated (cycling) myosin cross-bridges with actin, as well as, attached 
(non-cycling) latch-type cross-bridges; and (b) a model of force generation and mechanical coupling 
to the contractile filaments and their attachments to protein structures and the skeletal framework 
of the cell. The two subsystem models are tested independently and compared with data. Like- 
wise, the complete (combined) cell model responses to voltage pulse stimulation under isometric 
and isotonic conditions are calculated and compared with measured single cell length-force and 
force-velocity data obtained from literature. This integrated cell model provides biophysically based 
explanations of electrical, chemical and mechanical phenomena in cerebrovascular smooth muscle, 
and has considerable utility as an adjunct to laboratory research and experimental design. 

Keywords: Mechanistic model, intracellular calcium dynamics, myosin phosphorylation, active tension gen- 
eration 



I. INTRODUCTION 

Contraction in vascular smooth muscle cells can be ini- 
tiated by mechanical, electrical, and chemical stimuli. 
Several different signal transduction pathways activate 
the contractile mechanism, but all of which first lead to 
an increase in intracellular Ca^^ concentration [Ca^^J^. 
This increase in [Ca^^J^can arise from: (a) Ca^^ influx 
through voltage-dependent sarcolemmal Ca^^ channels, 
and/or (b) a release of Ca^^ from internal Ca^^ stores 
[e.g., from the sarcoplasmic reticulum (SR)]. In the cy- 
tosol, free Ca^"*" binds to a special Ca^^ binding pro- 
tein called calmodulin (CM), and the calcium-calmodulin 
complex (CaCM) activates myosin light chain kinase 
(MLCK), an enzyme that facilitates the phosphorylation 
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of myosin light chains (MLC) in the presence of ATP. 
This MLC phosphorylation leads to the cross-bridge for- 
mation and cycling of the myosin heads along the actin 
filaments, generating the active force necessary for mus- 
cle contraction. The stressed actin filaments are coupled 
via a viscoelastic system to the cell membrane, resulting 
in changes in cell length. Experimentally, the macro- 
scopic contractile properties of a single smooth muscle 
cell are usually measured in terms of length-force (L:F) 
and force-velocity (F:V) relationships. 

Several models have been developed to investigate the 
particular aspects of cerebrovascular smooth muscle func- 
tion, for example, including membrane electrophysiol- 
ogy [1] , cytosolic calcium regulation , phosphorylation 
of myosin [3, 01 or mechanical behavior [^). However, 
none has sought to integrate these components into a 
single functional model. In this study, we develop an in- 
tegrated model of vascular smooth muscle function based 
on the electrophysiological, biochemical and mechanical 
data available on isolated smooth muscle cells. 

Figure[T]shows a functional block diagram of the model 
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FIG. 1. Functional Block Diagram of the vascular smooth muscle cell. The model consists of two coupled subsystems 
(electrochemical and chemomechanical). The membrane model describes the electrical response of the cell membrane to 
stimulation and especially the L-type Ca^""" current (/ca,L), which brings Ca^^ into the cell. The intracellular fluid compartment 
model is based on material balances of Na^, and Ca^"^, whereas the SR model contains a Ca'^'*' -induced Ca'^'*' -release (CICR) 
mechanism that releases Ca^"*" into the cytosol, and actively recovers it from the cytosol. The activated calcium-calmodulin 
(CaCM) drives the phosphorylation of the contractile kinetics model and subsequently active force development by the cell 
mechanics model. AMp and AM represent phosphorylated and nonphosphorylated states of MLC that forms a crossbridge 
attachment to the actin filament (See the main text for details). 



for the vascular smooth muscle cell. The whole model is 
composed of two specific subsystem models. The first 
is an electrochemical model that combines a Hodgkin- 
Huxley (HH)-type model of membrane electrophysiology 
with models of the intracellular fiuid compartment and 
the sarcoplasmic reticulum compartment. The mem- 
brane model contains descriptions of the known ionic 
membrane currents active in cerebrovascular smooth 
muscle cells. Material balance equations for cytosolic 
ionic concentrations of Na^, K+, and Ca^"^ provide the 
basis of the fluid compartment model under the assump- 
tion that the concentrations of these ionic species in the 
extracellular medium are held at constant levels (see Ta- 
ble |T] for numerical values), and that Ca^^ concentration 
in the cytosol is buffered by the Ca^^ binding protein 
calmodulin (CM) and other Ca^^ binding proteins. The 
material balance for Ca^^ also considers the Ca^^ fluxes 
entering and leaving the cell cytosol via the SR and sar- 
colemmal membranes. The experimental data used to 
validate the electrophysiological properties of the mem- 
brane and fluid compartment models comes largely from 
voltage clamp experiments on isolated cells. This sub- 
system model thus characterizes the electrochemical be- 
havior of the cell and most importantly the regulation of 
intracellular Ca^^ . 

The second subsystem model is a chemomechanical 
model that has two coupled units: (i) a biochemical 
model of CaCM-dependent myosin light chain phospho- 
rylation and attached cross-bridge formation kinetics, 
and (ii) a mechanical model of force generation and me- 
chanical coupling within the cell. A modified version of 
the multiple-state kinetic model of myosin phosphoryla- 
tion originally developed by Hai and Murphy 3i, 4;] is used 
to represent the dynamics of myosin-actin cross-bridge 
attachment and detachment, and a mechanical model is 
used to describe the active force generation by the con- 
tractile filaments and the viscoelastic properties of their 
coupling and attachment to the cell wall. We use data 
from the length-tension studies of Fay and Warshaw @- 



|10[ on isolated smooth muscle cells to help validate this 
chemomechanical model. 

II. MODEL DEVELOPMENT 

A. Electrochemical model 

Figure [5]^ represents a lumped Hodgkin-Huxley type 
electrical equivalent circuit of the smooth cell mem- 
brane, which consists of a whole-cell membrane capac- 
itance (Cm) shunted by a variety of resistive transmem- 
brane channels, as well as, ionic pump and exchanger 
currents. The Kirchhoff's current law applied to the cir- 
cuit of Fig. yields the following differential equation 
describing changes in the transmembrane potential Vm- 

—77- = —7;— (-^Ca,L + Ik + ^K,Ca + ^Ki + 

Im + -^NaCa + ^NaK + -^CaP + Ib) (1) 

The individual membrane currents involved in this equa- 
tion are discussed in the subsequent sections. 

To account for temporal and voltage-dependent 
changes in the ionic concentrations of Na^, K+ and 
Ca^^ in the cytosol, a fiuid compartment model is also 
(Fig. 03) developed. Cytosolic material balance equa- 
tions (see Appendix I, Table HIl are given for Na+, Ca^+ 
and K"*" . Concentrations for these ions are assumed to be 
constant in the bulk extracellular bathing medium. The 
material balance equation for Ca^"*" accounts for Ca^"*"- 
binding to the cytosolic protein calmodulin and other 
nonspecific buffering media, Ca^^ uptake and release by 
the sarcoplasmic reticulum. 

1. Voltage-gated calcium current /ca,L 

Both L-type and T-type Ca^"*" currents have been de- 
scribed in whole-cell voltage-clamp studies on vascular 
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FIG. 2. The electrochemical Model. A. Membrane model 
describing ionic membrane currents and transmembrane po- 
tential. Black and white indicates resistance with voltage- 
dependent nonlinearity and the all white resistor indicates 
linear element. B. Fluid compartment model describing ionic 
dynamics, Ca^"*" buffering, and Ca^^ handling by sarcoplas- 
mic reticulum (SR). 



smooth muscle Regarding the time scales of 

Ca^^ channel activities, the L-type stands for a long- 
lasting inactivation, while the T-type stands for a tran- 
sient activation and inactivation. However, unlike the 
L-type calcium current, the T-type has been found in 
some arterial smooth muscle cells, but not all 17]. Single 
channel studies indicate that conductance of the L-type 
channel is approximately three times larger than that of 
the T-type channel (T-type, 8 pS; L-type, 25 pS [H]). 
Our study is focused on cerebrovascular smooth mus- 
cle cells and for these cells, we assume that the T-type 
channel is weakly expressed or is absent for the follow- 
ing reasons: (1) Voltage clamp studies of isolated vas- 
cular smooth muscle cells indicate that activation and 
inactivation of T-type Ca^^ current /ca,T occurs over a 
voltage range that lies 20-40 mV negative to the thresh- 
old for the L-type channel current /ca,T \)M j which sug- 
gests that /ca,T may not be responsible for sustained 
smooth muscle contraction. (2) The amplitude of T-type 
current has been reported to be around 10-20% of that 
of /ca,L and it is insensitive to dihydropyridine (PHP ) , 
which strongly inhibits maintained arterial tone |17| . (3) 
The T-type channels are found in relatively high den- 
sity in spontaneously active smooth muscle types and 
are associated with generating action potentials [IJ] , but 
in contrast action potentials are not generated in cere- 
brovascular smooth muscle cells, suggesting that the role 
by the T-type channels were minimal. 

The current /ca,L is assumed to be the only current 
underlying sustained Ca^"*" entry and maintained con- 
traction in cerebrovascular smooth muscle cells. Due to 
the strong voltage-dependence of /ca,L, a membrane de- 
polarization can modulate Ca^"*" entry and consequently 
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FIG. 3. Characteristics of the L-type calcium channel. A. 
The steady-state activation (di) and inactivation (/l) gating 
variables. B and C. Voltage-dependent time constants Td and 
Tfi for activation and fast inactivation, respectively. D. Rela- 
tionship between the peak 7ca l and the membrane potential. 
Data (*) from Rubart et al. [ll[| and (o) from Langton [l^ . 
E. Model-generated fits to measured voltage-clamp data from 
Rubart et al. [0| . Voltage clamp (for a duration of 200 ms) 
steps are from -30 to 30 mV in 10 mV steps. The holding 
potential is -60 mV. 



increase [Ca^^];. At potentials in the range -60 to -40 
mV, steady-state currents could be detected for min- 
utes [IIII13I, which suggests that /ca.L niay provide sus- 
tained Ca^"*" influx in the physiological range of mem- 
brane potentials above the resting state. 

Figure compares model-generated and measured 
steady-state voltage-dependent activation [Til, [11] and in- 
activation characteristics for the /ca.L- The model equa- 
tions used for /ca,L are given in Appendix I: Table IIIIl 
where the voltage- and time-dependent activation and 
inactivation variables are defined as and /i, respec- 
tively. Figs. [3J3 and C show the voltage-dependence of 
the associated activation and inactivation time constants, 
respectively. These time constants are modeled by Gaus- 
sian functions with an offset. We employ two time con- 
stants in describing the inactivation process for /ca,L, 
since Rubart et al. [ll| report that both a fast and a 
slow component of inactivation may exist, with the slow 
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time constant having a magnitude larger than 2 seconds 
or having a partial non-inactivation. Data obtained using 
200 ms voltage clamp pulses [ll| does not show long term 
current inactivation. Therefore, we assume that the /ca.L 
current has a non-activation (Table Fig. [3jl) shows 
peak instantaneous /ca,L as a function of membrane volt- 
age. This current is maximum at 10 mV and reverses its 
polarity at 60 mV. 

The parameters associated with the voltage-dependent 
steady-state activation and inactivation functions and 
the associated time constants are identified by us- 
ing a nonlinear least square optimization method (i.e., 
Levenberg-Marquardt algorithm) ^9i]. Fig. [3j5 compares 
model-generated and experimentally measured voltage- 
clamp dataset of Rubart et al. [ll| for /ca,L- 



2. Stretch-sensitive current Im 

Stretch-sensitive ionic membrane currents have been 
shown to be present in arterial smooth muscle cells [ij, 
[TgI . [2oI | . When activated, these currents produce signifi- 
cant changes in membrane potential. In vascular smooth 
muscle cells, the stretch-sensitive membrane current Im 
has a reversal potential about -15 mV [20]. We assume 
that the lumped channel conducting Im has ionic com- 
ponents of K+, Na^ and Ca^'^ (with selectivity order: 
K+ > Na^ > Ca^'^), and we express each of these com- 
ponents using a modified Goldman-Huxley-Katz equa- 
tion: 



Im — Im,K + ^M,Na + ^M,Ca 



^M.K 



lM,Na — PP'Nz 



[Na+]o - [Na+],e^-^/«^ 

1 _ f,V^F/RT 



, [Ca^+]o - [Ca^+].e^-W 

iM,Ca — PZ fCa. 



1 _ ^V^zF/RT 



1 ^ g-(0"-0"l/2)/'^ 



(2) 
(3) 

(4) 

(5) 
(6) 



where the coeflticient /3 = AmPmF'^Vm / {RT) . The quan- 
tity Am is the surface area of the smooth muscle cell. 
We further assume that: (1) the geometry of the cere- 
brovascular smooth muscle cell is cylindrical; (2) the cell 
volume is 1.0 pi; and (3) the length of a cell is 150 //m. 
These geometric values can be used to calculate the cell- 
surface area Am- The constants Pk, P^a. and Pca rep- 
resent the permeability of the lumped channel to K+, 
Na^ and Ca^^, respectively. The permeability ratios for 
K+, Na+ and Ca^+(3.0:2.0:1.0) are adopted from Davis 
et al. |20j, and the open channel probability P,„ is char- 
acterized by a Boltzmann's relationship that is a function 
of the stress a experienced by smooth muscle membrane 
(ci/2 = 204 mmHg and fco- — 2.03 mmHg). The stress a 
is measured as force per unit cross-section area. Davis et 
al. reported that a 10-15% stretch leads to a 40-50 pA 




FIG. 4. Characteristics of voltage-dependent potassium chan- 
nel (i.e., the delayed rectifier). A. Steady-state activation gat- 
ing variable pK c omp ared to data from (*) Volk et al. [2ll] and 
(o) Nelson et al. [23| B. and C. Activation time constants Tpi 
and rp2, respectively, compared with data (*) from Volk et 
al. [211 ] D. Model generated fits to measured voltage-clamp 
data (Volk et al. [2l|].) Voltage clamps (with a duration of 
400 ms) are applied from the holding potential (-60 mv) to 
10-60 mV in steps of 10 mv. 



whole-cell inward current through stretch-sensitive ion 
channels 12011. 



3. Ca.^ -activated K+ current 

In vascular smooth muscle cells, the outward current 
is attributed lar gely to the Ca^'^-activated K+ current 
(-^K.Ca) [H, m, Hij. This current is activated in re- 
sponse to a large influx of Ca^'''as the result of mem- 
brane depolarization. An increase in /K,Ca brings about 
a hyperpolarization due to K"*" efflux, which limits the 
voltage-dependent Ca^"*" , and leads to smooth muscle 
relaxation. Patch clamp studies [l^l have shown that 
single Ca^'^-activated K"*" channels in smooth muscle cell 
isolated from rat cerebral arteries have a mean conduc- 
tance of 207 pS. The time course of channel opening con- 
sists of two distinct components, i.e., fast and slow Ca^^- 
activated processes with mean time constants of 0.5 and 
11.5 ms, respectively. The steady-state open probabil- 
ity of the channel is well-described by Boltzmann rela- 
tionship that is shifted leftward by a 45mV per decade 
increase in [Ca^''']i(see Table Hvl for its mathematical de- 
scription). 



4. Inward rectifier current Jk; 

The inward rectifler current (/ki) plays an important 
role in establishing the restin g m embrane potential of ar- 
terial smooth muscle cell psi. |25|| . Control of channel ac- 
tivity by the extracellular potassium concentration [K+]o 
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FIG. 5. A four-state kinetic model of ryanodine receptor ac- 
tivation controlling Ca^^-induced Ca^^ release (CICR).The 
receptor has two regulatory sites that bind to Ca'^'*' for acti- 
vation and inactivation. 



constant associated with inactivation process is relatively 
large (2910 ms at K„=60 mV, Volk et al. [2l|), we do not 
include this slow process in our model. The activation 
process is characterized by two exponential components 
with different time constants. The voltage-dependent 
steady-state activation function is shown in Fig. 
whereas the time constants of activation Tpi and Tp2 are 
shown in Figs. |3J3 and C, respectively. Model-generated 
fits to measured voltage clamp data (Volk et al. [2l|) are 
shown in Fig. Up. The half-activation voltage (Vi/2,/f) 
and slope factor (k) of the Boltzman function, as well 
as, the parameters associated with voltage-dependent 
time constants are determined by using a non-linear least 
square parameter estimation method [l^. The values 
used in the model for Vi/2./f (-1.77 mV) and k (14.52 
mV) are in good agreement with measurements by Volk 
et al. f2l| {Vi/2,K = -6 mV, k = 9 mV) and by Nelson 
et al. (^1/2,^ = -9 niV, A: = 11 mV). 



is also a distinguishing feature of this current. Studies by 
Quayle et al. [26| characterize the maximum slope conduc- 
tance as a function of [K~^]o according to: 

5max,Ki -Gk.([K + ]o)'*' , (7) 

where Gki and are constants with typical values 
of 0.145 and 0.5, respectively for pig coronary arterial 
smooth muscle [26|. We have used a similar relationship 
in our model of rat cerebral arteries, and have adjusted 
the parameters of this equation to better characterize 
measured data from rat cerebral arteries [25j . The pa- 
rameter values used are given in Table HI The chord con- 
ductance gK — I /{V — Ek) has a Boltzmann relationship 
with membrane potential, and the mid-potential Vi/2,k. 
is a function of [K+]o . To characterize this relation- 
ship, we used the mathematical representation proposed 
by Quayle et al. Q 

Vy2M.^A\og^^[K% + B . (8) 

Typical values for A and B are 25.19 mV/log]^o(™-'^) 
and 112.29 mV, respectively. When the extracellular 
concentration [K"'']o increases, the Boltzmann curve is 
shifted rightward toward depolarized potentials (see Ap- 
pendix IHl Table |Vl|. 

5. Delayed rectifier current Ik 

The delayed rectifier current (Ik) activates when the 
membrane is depolarized, and most vascular smooth mus- 
cle cells have this current. It provides an outward current 
in response to membrane depolarization, which is impor- 
tant to counterbalance the effect of inward currents, thus 
limiting further membrane depolarization. 

In our model, the steady-state activation process of 
Ik is voltage-dependent, and it is described by a Boltz- 
mann function (see Appendix I: Table IV)) . Since the time 



6. Na^/Ca^^ exchanger and membrane pumps 

Our model of the sarcolemma includes membrane cur- 
rents associated with the Na"''/Ca^~'' exchanger (/NaCa), 
the Na'''/K+ pump {lysaK ), and the Ca^'''pump (/cap)- 
Na~''/Ca^''' exchange is a counter-transport system that 
translocates Na"*" and Ca^^ across the plasma membrane. 
This transporter has been well investigated and modeled 
for cardiac muscle cells [2^. Studies of this system in 
smooth muscle cells have also been reported (28|. The 
Na~''/Ca^''' exchange system mediates either Ca^"*' efflux 
in exchange for Na"*" uptake or Na"*" efflux in exchange 
for Ca^''" uptake depending on membrane potential (Vm) 
and the gradients for Ca^"*" and Na"*". The stoichiometry 
of the exchange appears to be 3Na+ : lCa^+ El. 

As an active transport system, the Na^/K+ pump is 
stimulated by the Na~'~/K+ ATPase and powered by hy- 
drolysis of ATP to move Na"*" and K+ ions against their 
electrochemical gradients. The stoichiometry of this ac- 
tive transport system is SNa"*" : 2K"'" and its ion transfer 
characteristic has a sigmoidal dependence on membrane 
potential [H. The Na+/K+ pump is electrogenic and 
provides a net outward current that helps to maintain 
the Na^ and K+ concentration gradient across the cell 
membrane. 

Activity of the Ca'^^ pump is regulated by calmod- 
ulin (CM), which serves to stimulate Ca^"*" extrusion and 
Ca^+-ATPase activity ^28^. In our study, we use a simple 
sigmoidal relationship to model the coupling between the 
whole cell Ca^^ pump current and intracellular concen- 
tration of cytosolic Ca^^, and assume that the CM gat- 
ing kinetics of the Ca^^ pump are very fast so that the 
binding process can be considered instantaneous (see Ap- 
pendix: Table rvTI)) . Some studies report that the /NaCa 
does not play a significant role as a Ca^^-removal mech- 
anism in the smooth muscle cells of rat cerebral arter- 
ies [i3,H[3H. 
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7. CICR mechanism of sarcoplasmic reticulum 



In cerebrovascular smooth muscle cells, there is no 
general agreement regarding the importance of the sar- 
coplasmic reticulum (SR) in regulating [Ca^''']i. One 
group maintains that the Ca^"*" transient is the re- 
sult of sarcolemmal /ca,L influx, and suggest that lo- 
cal Ca^'''-release from ryanodinc-sensitive Cs?'^ into con- 
fined spaces regulates membrane potential via activation 
of -^K.Ca [H, [s^. Others maintain that, in addition. 



there is a significant contribution to global intracellu- 
lar [Ca^+], by SR Ca^+-release [sl. We model Cs?+- 
induced Ca^^-release (CICR) using a multiple-state ki- 
netic model, which is similar to models developed for the 
SR of the cardiac ventricular cell (35l - [37| . As shown in 
Fig. O the gating of a lumped ryanodine-sensitive Ca^^ 
release channel is regulated by the Ca^~''-dependent re- 
ceptor kinetics. We assume that the secondary release of 
Ca^^ by SR affects the cytosolic Ca^"*" concentration di- 
rectly. The following differential equations describe this 
dynamic behavior of the different states of the ryanodine 
receptor (RyR): 



dK 

dt 
dRm 
^ dt 
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(9) 



subject to the mass conservation constraint: 



are expressed as: 



R 



00 



i?oi + -Rio + -Rii — 1 , 



(10) 



where, the four different states of ryanodine receptor frac- 
tions are denoted as Rqq: free receptors; Riq: recep- 
tors with Ca^"*" bound to activation sites; Rqi- receptors 
with Ca^^ bound to inactivation sites; and Rn: recep- 
tors with both activation and inactivation sites bound by 
Ca2+. The open probability of the Ca^"*" release channel 
controlled by this receptor mechanism is governed by the 
state variable Riq. The rates of Ca^"*" binding to activa- 
tion and inactivation sites are independent. Binding to 
activation sites is much faster than to inactivation sites 
(e.g., i^ri is much larger than Xr2). In the model, there 
are two Ca^'^ sites used for activation and one for inacti- 
vation and a second-order gating of release currents (/rd) 
by Rio is applied. The above four ordinary differential 
equations (jO)) finally reduce to a 3rd order system due to 
the algebraic constraint given in Eq. pH)) . 



The SR itself is modeled as two-compartment (re- 
lease/uptake) store with a fluid connection. A diffu- 
sion current Itr flows between each the compartments, 
as shown in Fig. [5) The SR is refllled via active pump- 
ing of Ca^"*" into an uptake compartment, whereas the 
release compartment is responsible for storage of Ca^^ 
prior to release into the cytosolic space upon activation 
of the RyR. Intra-store Ca^^ balances and current fluxes 



[Ca- 



2+1 



[Ca^+l 
([Ca2+]„ - 



m,up 
2 



[Ca^+],)(2^^-vol„) 



Ttr 



(11) 

(12) 



d[Ca^ 



Irel — RfoiiCB, 
]up -^up -^ti 



2+1 



ic.-,.,e^3^(i3) 



Trcl 



dt 



dt 



2F ■ vol„ 
2F ■ voir 



(14) 
(15) 

where the uptake current /up is regulated by Ca^"*"- 
ATPase activity, and the Ca^"*" release current I^ei is gated 
by the activation state i?io of ryanodine receptor model 
and the Ca^"*" concentration gradient across the SR mem- 
brane. Description and numerical value of parameters 
associated with the model are listed in Table IIXI 



8. Ionic balances of intracellular materials 

We have also developed a fluid compartment to ac- 
count for changes in the cytosolic concentrations of Na"*", 
if"*" and Ca^"*". This model includes a description of the 
geometry of smooth muscle cell, as well as Ca^'^ buffer- 
ing by calmodulin (CM). Since the calcium-calmodulin 
complex (CaCM) plays such a vital role in regulating 
the myosin phosphorylation process, we distinguish this 
buffering process by modeling it separately. Ca^^ buffer- 
ing by other elements such as the mitochondria, and the 
fluorescent indicator dye sensitive to [Ca^"'']i, are lumped 
together as "other" buffering processes ([i?F], equation 
P^ ). Material balance equations for intracellular ions 
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FIG. 6. A four-state kinetic model of CaCM dependent 
myosin phosphorylation and cross-bridge formation. This fig- 
ure is adopted from Hai and Murphy [3,0]. 



are given in Appendix: Table |TI1 The extracellular ionic 
concentrations are assumed to be constant ([Na^]o = 140 
niM, [K+]o = 5 niM and [Ca^+]o = 2 niM). Calmod- 
ulin has four free sites for binding Ca^"*". The average 
binding affinity of calmodulin sites for Ca^"*" can be ex- 
pressed as a dissociation constant (Kd) of 2.6x10""^ M 
based on the assumption that there are two high-affinity 
sites with Kd =2.0xl0~^ M, and two low-affinity sites 
with Kd =2.0x 10-6 M [H, m. The Ca^+ buffering 
by calmodulin and other Ca^^ binding {[Bp]) elements 
in the cytosol are described as first-order dynamic pro- 
cesses, which are represented by the following differential 
equation: 



d[ScM] 
di 

[CaScM] 

dt 
[CaBF] 



k^iiCaScu] - ki[Ca?+UScM] (16) 

[5cm] - [5cm] (17) 

k^d[CaBF]^kd[Ca''+UBF] (18) 

[Bp] - [Bp] , (19) 



where [^cm] represents the concentration of free calmod- 
ulin sites for Ca^"*" binding and [S'cm] represents the total 
concentration of calcium binding sites available for Ca^"*" . 
The constants k-i and ki are kinetic rate constants for 
Ca^"*" uptake and release, respectively, and the ratio of 
k-i to ki is equal to dissociation constant. Similar anal- 
ysis is applied to the Ca^"*" buffering process by other 
nonspecific proteins as described in Eq. p^. 



B. Kinetic model for smooth muscle cell 
contraction 

Interaction between actin and myosin filaments pro- 
vides the molecular basis for smooth muscle contraction 
which is regulated by intracellular Ca^^. Contractile 
mechanism of smooth muscle exhibits features which are 
similar to as well as distinguishing from striated mus- 
cle. A strong difference in smooth and striated muscle 
is that tension can be maintained in smooth muscle cells 
when [Ca^~'']iand the level of phosphorylation are rela- 
tively low. This observation has led to the latch-bridge 
hypothesis of vascular smooth muscle contraction |4C| . 
Although other mechanisms may co-exist for the regula- 
tion of contractile strength in smooth muscle {e.g., pro- 
tein kinase C may activate force generation with small or 
no changes in [Ca^''"]i[41|]), we assume that CaCM-MLCK 
dependent activation is the only mechanism present, and 
employ the four-state kinetic model developed by Hai and 
Murphy (H-M) [3, 4] to describe myosin phosphorylation 
and latch bridge formation. 

Figure [6] is a depiction of H-M model, which consists 
of four fractional species: free cross-bridges (M), phos- 
phorylated cross-bridges (Mp), attached phosphorylated, 
cycling cross-bridges (AMp) and attached dephosphory- 
lated non-cycling cross-bridges (as latch bridges (AM)). 
The latch state has been found to be unique in smooth 
muscle, which endows the smooth muscle cell with the 
ability to maintain force when the level of calcium- 
dependent myosin phosphorylation is relatively low. The 
kinetics of these four species can be described by the fol- 
lowing differential equations: 



/ dU/dt \ 
dMp/dt I 
dAMp/dt 
V dAM/dt J 





~Ki 


K2 







/ M 


Ki 


-K2 - K^ 


K4 


» ] 


Mp 







K3 


-K4 - 




AMp 












-Ke - K7J 


\AM 



(20) 



subject to the constraint 

M Mp + AMp + AM = 1 , (21) 

where the fraction of phosphorylated myosin is defined as 
the sum Mp-f AMp, whereas attached cross-bridges are 
represented by the sum AMp+AM. In our modification 



of H-M model, the CaCM-dependence of myosin phos- 
phorylation, rather than direct dependence on Ca^"^, is 
represented through rate constants Ki and Kq with an 
analytic relationship as described in Eq. (22). Ki = Kq 
is assumed, which means that the phosphorylation pro- 
cess has same rate in changing from M to Mp, as from 
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AM to AMp. Thus, 



[CaCM]' 



[CaCM]2 



(22) 



where iiTcaCM is the half-activation constant for the 
CaCM-dependent phosphorylation rate constant. This 
sigmoidal relationship indicates the saturation of the rate 
constants Ki and Kq when the concentration of CaCM 
is large. We also assume the dephosphorylation process 
from Mp to M occurs at the same rate as from AMp to 
AM, I.e., K2 = K^. 



C. Mechanical properties of the cell 

Experimental studies conducted by Warshaw and Fay 
on the length-force relationship of the single smooth mus- 
cle cell report that cross-bridge stiffness is between 1.2 
and 1.5 times greater than the measured cell stiffness. 
They report that the elastic response originating within 
the muscle cell is attributed to cross-bridges acting in 
series with an elastic element that has an exponential 
length- force relationship [1, . A viscous contribution to 
cell stiffness has also been detected in experiments where 
the phase shift between an applied cell length change and 
the observed change in force has been measured [8, 9]. 

A modified Hill model [42[ is used to describe the me- 
chanics of the cell (Fig. [7]). This model reflects the cou- 
pling between active the force generation by the cross- 
bridge mechanism and the mechanical coupling proper- 
ties of the myofilaments and their viscoelastic attachment 
to structures in communication with the cell wall. 



1. Passive Force 

The passive elasticity of the cell {kp) is modeled by an 
exponential force-length relationship used to represent 
the nonlinearity of the passive elasticity. The passive 
force is given by: 



1 



(23) 



where kp is a constant and Iq is the length at which the 
cell generates zero passive tension. 



2. Cross-bridge elasticity 

Due to the elasticity of attached cross bridges, cell 
tension can also be sustained by cross-bridges, and this 
length-dependent force is characterized by: 

= {k^, ■ AMp + k^, ■ AM) . • e-'3('"/'-'-i)' , (24) 

where is the extension of cross-bridge and la is the 
length of active element which indicates the overlap be- 
tween the myosin and actin filaments. The optimal 



length (/opt) of the active element is the length at which 
maximum overlap occurs between the myosin and actin 
filaments, and it is also the length at which the cell de- 
velops maximum tension for the active state achieved. 
We approximate the effect of filament overlap in the de- 
velopment of active tension by the Gaussian function in 
Eg. d^^ . This equation for also reflects that both 
types of cross-bridges (cycling and latch bridges) con- 
tribute to the elasticity of the attached cross-bridges. 
The spring constants k^i and fcx2 denote the maximal 
stiffness that can be achieved by phosphorylated at- 
tached cross-bridges and latched cross-bridges, respec- 
tively. States AMp and AM are the outputs of the kinetic 
contractile model described in previous section, and in- 
dicate the distribution of attached cross-bridges and the 
state of interaction between myosin and actin filaments. 
Consequently, the effective stiffness contributed by phos- 
phorylated cross-bridge is the product fcxi AMp, whereas 
that contributed by latch bridges is fcx2 AM. Eq. ([M)) 
also indicates that the two types of attached cross-bridges 
with different effective stiffness contribute in parallel to 
the total active force and the total stiffness is further 
modulated by cell contractile kinetics which result in a 
change in active element length la. 



3. Active force generation 

The active force generated by attached cross-bridge is 
given by: 



Fa^ 



fAMp • AMp(i;, + ^) 



fAM-AM^ 

at 



-I3{la/la 



(25) 

where /amp and /am are friction constants for phospho- 
rylated cross bridges and latch bridges, respectively. In 
this formulation, we assume: (1) Only attached phos- 
phorylated cross bridges exhibit cycling behavior which 
produces translational motion of the actin filament at a 
velocity v^, and this is the only process contributing to 
active force generation. Latch bridges lose cycling capa- 
bility, and their continuing attachment to the actin fila- 
ment increases the stiffness of the cross-bridge coupling 
[Eq. (p4)) ]: (2) Both species of attached cross bridges, 
phosphorylated and latch bridges, slide along the actin 
filament with velocity dla/dt in response to cell stretch, 
release and contraction. 

The first assumption indicates that populations of at- 
tached cross-bridges are heterogeneous, and that latch 
bridges act as an internal load on fast cycling cross- 
bridges m, nil . Note that both the total effective stiff- 
ness and the friction offered by the cross-bridges are sub- 
ject to modulation by la, the length of the active element, 
through a Gaussian relationship that is expressed in Eqs. 
([M)) and (|25p . respectively. The Gaussian function ex- 
presses the effect of the degree of overlap of the actin 
and myosin filaments. 
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FIG. 7. Mechanical model for the smooth muscle cell. The passive length-force properties of the cell are described by the 
spring with stiffness kp, whereas the springs characterized by stiffness k^c and ks denote the elasticity of attached cross bridge 
and series elastic element, respectively. The wall attachment coupling also has a viscous element characterized by the damping 
coefficient Hs- The active contractile element (length la) models cross-bridge actin filament interaction by the forces generated 
by cross-bridge cycling and sliding. 



4- Series viscoelastic force 



III. RESULTS 



The series coupling element is modeled as an elastic 
Kelvin- Voigt body that consists of a spring (by an elas- 
tic constant ks) connected in parallel with a dashpot (by 
a frictional constant Hs)- The total force across the body 
is given as the sum of elastic and viscous forces, where 
the spring is characterized by an exponential length-force 
relationship, and the force sustained by the series vis- 
coelastic element is expressed as: 



dl 



(26) 



where the viscous coefficient /is and kg are constants, and 
Iso is the length at which the series elastic component 
sustains zero force. 

At each instant of time, the force sustained by (a) the 
cross bridges, (b) the active force generation element, and 
(c) the series element is the same, i.e., 



F — F — F 
The total force sustained by the cell is: 

Ft^Fa+Fp . 



(27) 



(28) 



The numerical values of parameters associated with this 
cell mechanics model are listed in Table [X] 

Equations and parameters of the electrochemical 
model are listed in Table U through Table IXll Initial con- 
ditions for typical model simulations are listed in Ta- 
ble IXIII The complete model consists of 23 state vari- 
ables. Model parameters are identified by using numer- 
ical optimization methods with non-linear least square 
techniques [l^ to fit experimental data. The model is 
implemented and the results are visualized by using the 
MATLAB programming language on an Intel-based PC 
platform. 



We test the cell model by subjecting it to condi- 
tions that emulate various experimental protocols used 
in testing the functional behavior of single smooth mus- 
cle cells, and compare model predictions with corre- 
sponding measured data from the literature. Agreement 
between model predictions and measured data over a 
wide range of testing protocols indicates the capabil- 
ity of this model to represent key mechanisms involved 
in excitation-contraction coupling, as well as, the gen- 
eral electromechanical behavior of this cerebrovascular 
smooth muscle cell. 



A. Testing the electrochemical model 

1. Calcium transient 

The Ca^^ transient is of importance for proper smooth 
muscle function, since it represents the input to the con- 
tractile mechanism. In cerebrovascular smooth muscle, 
Ca^^ entry into the cytosol is mainly via the inward 
voltage-dependent membrane current /ca,L [13 ■ The 
principal Ca^^-removal mechanisms at work in this cell 
type are the separate and distinct sarcolemmal and SR 
Ca2+-ATPase pumps. Kamishima et al. [sij report that 
in the rat cerebral arterial smooth muscle cell, the rate of 
Gs?^ removal exhibited by cells bathed in a Na"*" free so- 
lution was not significantly different from that observed 
in control cells. Studies by Ganitkevich et a/.fs^ on 
smooth muscle cells of guinea pig urinary bladder, also 
indicate that Ca^"*" extrusion via /NaCa is negligible in this 
cell type. Although our model contains a mathematical 
expression for /NaCa, the magnitude of the sarcolemmal 
pump current /ca,p and the SR uptake current /up are 
made relatively large compared with /NaCa7 and hence 
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FIG. 9. CaCM-dependent myosin phosphorylation and cross- 
bridge distribution. Calcium transient is eUcited with 1.6 s 
voltage depolarization from -60 to mV. A. Waveform of the 
concentration of CM sites occupied by Ca^'*' . B. Total fraction 
of attached cross-bridge population ( AM-f AMp) . C. Tempo- 
ral distribution of attached cross-bridges in phosphorylated 
(AMp : solid line) and dephosphorylated states (AM: dashed 
line). 



FIG. 8. Calcium transient in response to a voltage pulse. 
A. Intracellular Ca'^^ transient in response to 1.6 s step de- 
polarization to mV from a holding potential of -60 mV. 
Model-generated data is plotted with measured data from 
Kamishima et al. [3ll . H^l). Here the change in [Ca^+]ifor 
model-generated and measured data are each normalized by 
the maximum peak value [Ca'^''']niax encountered in the series 
of voltage clamp tests (see Fig[Sl), which happens at 10 mV 
as shown in panel E. B. Voltage-dependent L-type calcium 
current /ca,L in response to voltage step (note, this panel has 
a different time scale than others). C. Calcium removal cur- 
rents: the sarcolemma Ca^'*'-ATPase pump current Ica.p and 
the SR membrane Ca^"'"-ATPase uptake current /up. D. Ca^"*" 
concentrations in the uptake and release compartments of the 
SR. E. Voltage dependence of change in [Ca^^]i(A[Ca^^]i) 
measured relative to baseline [Ca'^''']i. 



-^NaCa is Considered small. 

Calcium transients can be induced in several ways in 
smooth muscle cell. Here we focus on the voltage clamp 
method and emulate the experimental protocols used by 
Ganitkevich et al. [lOl and Kamishima et al. [3l|, [s^l to 
produce Ca^"*" transients. Fig. [8]A_ compares the model- 
generated Ca^"*" transient in response to a 1.6 s voltage 
pulse from holding potential of -60 mV to a depolarized 
level of mV. The measured data shown in this figure 
is from Kamishima et al. [Sll, • Figure [5]4. shows that 
the depolarization pulse triggers an increase in cytosolic 
calcium concentration ([Ca^^];), whereas Figs. [SJB and 
C show that the leading edge of this Ca^"*" transient is 
due mainly to /ca.L, but also to the SR release current 
/loi. Decay in the Ca^^'^transient is mainly caused by 
the rapid decline in /ca,L coupled with the increase in 
sarcolemmal and SR uptake pump currents (/caP and 



Jup, respectively). Fig. [8jl) shows that [Ca^^] within the 
SR compartments is quickly depleted with the onset of 
the voltage pulse and the SR Ca^+ release (/rci)- The 
model predicts that the SR compartments are slow to 
refill after the depletion of the store (Fig|8p). 

Experiments by Ganitkevich et al. [30| on single 
smooth muscle cells of the guinea-pig urinary bladder 
show that the voltage-dependence of peak Ca^"*" tran- 
sients resembles the bell-shaped peak I-V relationship of 
Ica,L (Fig-[5P)- Fig. [5J5 shows that the model-generated 
peak Ca^"*" transients exhibit a bell-shape relationship 
with membrane potential, similar to that exhibited by 
^Ca,L- Both curves have maxima at 10 mV. The protocol 
followed by Ganitkevich et al. [s^l to study the voltage- 
dependence of the Ca^"*" transient, involves the delivery 
of a number of 1.6 s depolarization pulses of different 
amplitudes over the range -40 mV < V < +60 mV in 10 
mV steps, starting at the holding potential of -60 mV. 
CytosoHc /ca,L and the [Ca^^Jitransient are measured at 
each voltage step. 

B. Testing the chemomechanical model 

1. CaCM-dependent myosin phosphorylation 

Increases of intracellular Ca^^ initiate the 
smooth muscle contraction through the postulated 
Ca^^— >CaCM— >MLCK pathway of myosin phospho- 
rylation and cross-bridge formation. Model simulation 
of this cellular event is shown in Fig. [HI A 1.6s mem- 
brane depolarization from -60 to mV is applied to 
generate the cytosolic Ca^"*" transient. Fig.lHK shows the 
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FIG. 10. Length-force relationship in the isolated smooth 
muscle cell. Model-generated data compared with measured 
data (*) is from Harris et al. [2]. 



dynamics of the concentration of the CM sites occupied 
by Ca^"*" ion in response to the [Ca^~'']itransient. Kinetic 
model simulation of CaCM activated myosin phosphory- 
lation give the response of total attached cross-bridges 
population as shown in Fig.lHB: whereas the distribution 
and transition of cross-bridges between phosphorylated 
state and latch state is shown in Fig. [9l3. The majority 
of attached cross-bridges at the resting level of -60mV 
are latch bridges. Phosphorylation of free myosin M and 
the transition AM— >-AMp following the Ca^''' transient 
contributes to the rapid increase of the active cycling 
cross-bridge population. As the CaCM level decreases, 
the dephosphorylation process becomes dominant and 
the attached cross-bridge population tends to transfer 
back into the latch state. 



2. Length-force relationship 

Classic isometric test protocols have been applied to 
study the static length-force (L:F) relationship in isolated 
smooth muscle cells Specifically, our model simula- 
tions employ a very slow ramp of length increase from 
Iq (see Appendix I VII Table |X| and the change in force 
is observed. In Fig. [TUl the L:F relationships predicted 
by the; jnodel are compared with measured data (Harris 
et al. [6] on toad stomach smooth muscle cells). For cell 
lengths below the optimal length (maximum overlap be- 
tween myosin and actin filaments; Lceii in Fig. 1101) . there 
is an approximately linear relationship between length 
and force. Figure fTOl shows that force tends to approach 
zero when cell length becomes less than 0.4 of Lceii, which 
agrees well with the measurements of Harris et al. Q . 



3. Force-velocity relationship 

A classic isotonic quick-release technique is frequently 
applied to the study of the relationship between force and 
shortening velocity in muscle. Measurements by War- 
shaw ^] on isolated single smooth muscle cells of toad 
stomach muscularis shows that force- velocity relationship 
of the single cell agrees with the well-known hyperbolic 
description of Hill equation [iH] : 

+ + h) = {l + , (29) 

-^max -^max -^max 

where Fmax is maximum isometric active force, and a and 
h are constants. 

In isotonic tests, the cell is electrically activated [7|. 
The output of the contractile kinetic model is the sum of 
the attached cross-bridges (AM -I- AMp). Compared to 
the duration of AM-|-AMp in response to a calcium tran- 
sient (typical duration: 50-70 s), the duration of the a iso- 
tonic test is much smaller (1.5-2 s). Thus, myosin phos- 
phorylation kinetics are assumed to be relatively con- 
stant during isotonic conditions. Consequently, AM and 
AMp are set to constants during isotonic test procedures 
(typical values AM==0.3 and AMp=:=0.5; Fig.[l]). Model- 
generated length and velocity responses to quick releases 
in force are shown in Fig. 1111 The maximum isometric 
force Fjnax is the force sustained by the cell before the 
initial force release, and it is a function of the contrac- 
tile kinetics and cell length (Lccii). -Fmax is used as the 
reference to a series of force releases, which have various 
magnitudes ranging from 0.1 to 0.9 i^max with a step-size 
0.2 fmax- Force release to the test level is complete within 
50 ms for each trial. The test force level is maintained 
for 500 ms before the cell is brought back to i^max over 
a 1.0 second period (Fig. [TTR). Cell length and velocity 
responses are shown Figs. [TTB and C, respectively. 

Figure HID sdiows model-predicted and measured data 
by Warshaw Q for the force-velocity (F:V) relation- 
ship of the single smooth muscle cell. Warshaw reports 
that the F:V relationship can be fitted by Hill hyper- 
bola as described in Eq. (f29| with a/Fj^ax =0.268 and 
b = 0.163Lceii/s Ji]. Panel D also shows that our model 
can produce reasonable fits to measured F:V data, in that 
it provides good fit to the Hill equation, which was used 
by Warshaw [3j to fit their data. 



C. Complete model testing 

The full or complete model is now subjected to tests 
that highlight the relationships between membrane po- 
tential, cell strain, calcium transients force generation 
and mechanical response. Two stimulation protocols are 
applied: (1) membrane depolarization induced by voltage 
pulse and (2) strain change induced by cell stretch. 
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Lceii is constant (123 fim). Importantly, Figs. [T2F-H 
show that peak Ca^"*" corresponds directly to the peaks of 
total force (F) and the rate of force development {dF/ dt) , 
which is consistent with measured data [13, . 

In addition, we investigated the whole cell response un- 
der isotonic conditions. In this case, the modeled smooth 
muscle cell is subjected to the same electrical stimulation 
protocol delivered under isometric conditions, but rather 
than a length clamp, a constant force {F = 0.4 /LtN) is ap- 
plied to the cell. The electrically elicited Ca^~''transient 
is shown in Fig.fT^K. which also resembles the Ca^"^ tran- 
sient shown in Fig. [T2K. Active contraction follows the 
Ca?^ transient and the change in cell length is shown in 
Fig. n^B. Contraction originates in the active force ele- 
ment of the mechanical model, and as the result of the 
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FIG. 11. Length and velocity responses in isotonic tests of 
the model and force-velocity relationship. A. Isotonic force 
release protocol: release magnitudes from 0.1 to 0.9-Fmax with 
step-size 0.2_F'max. Each force level is maintained for 500 ms. 
B. and C. Length and velocity changes in response to the 
force changes, respectively. Note that shortening velocity in 
F:V relationship should be the velocity taken after the iso- 
tonic release. -Fmax = 2.13 fiN, and Icdi ~ 123 fim. D. Solid 
line: model simulation of F:V relationship. Data (o) from 
the Hill hyperbola [see Eq. (|29p ] which was used by War- 
shaw to fit data from on toad stomach muscularis with 
a/Fmax =0.268 and b = 0.163 Lcdi/s. 



1. Electrical forcing 

Figure [T^ shows a complete model simulation with the 
cell held under isometric conditions. A depolarizing 1.6s 
voltage pulse from a holding potential of -60 mV to 
mV is applied, and the active force generated by smooth 
muscle contraction follows [Ca^"*"] transient with a delay 
fFig.fl^K and B). The voltage pulse is of very short dura- 
tion relative to the duration of the Ca^^ transient. The 
model predictions agree well with the measured data of 
Yagi et al. [10] on single smooth muscle cells. Due to the 
active cycling of the phosphorylated cross bridges, the 
length of active element shortens (Fig. [T2jl)), and corre- 
spondingly the lengths of the other two elements {Is and 
Ix) increase (Panels C and E) to accommodate the in- 
crease in force. Under isometric conditions, cell length 
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FIG. 12. Complete model testing with voltage pulses under 
isometric condition. A. Ca^^ transient in response to 1.6 s 
voltage pulse from the holding potential -60 mV to mv. B. 
The active force developed in the cell in response to the Ca^^ 
transient. C, D and E. Internal length adjustment of the 
smooth muscle cell: Is (series elastic element), la (active ele- 
ment) and Ix (cross-bridge element), respectively. Zccii ~ 123 
fim. F. Peak [Ca^^JiUnder different levels of voltage pulses 
(from -40 to 60 mV with 10 mV step). G and H. Peak Force 
and peak rate of force development under different levels of 
voltage pulses, respectively. The fastest rate of force develop- 
ment corresponds to the largest peak value of [Ca^+]i. 
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cross-bridge cycling interaction with actin filament, the 
length of the active element la decreases to make contri- 
bution to the cell contraction (Fig. [l"3t)) At the leading 
edge of this response, the length of the series element and 
the extension of cross-bridge (Ig and Ix, respectively) in- 
crease, and because of the imposed isotonic condition, 
the force in this branch increases to counterbalance the 
force reduction in the passive branch due to the cell con- 
traction (Figs. fT3l C and E). Afterward, the cell begins 
to relax when the Ca^^ removal mechanisms (/caPi^up) 
takes over the cytosolic Ca^"*" regulation. 



2. Strain forcing 

This simulation tests the effect of a strain applied to 
the smooth muscle cell on the contractile mechanism. In- 
ward current through the stretch-sensitive channel (Im) 
plays a critical role in depolarizing the membrane, which 
regulates the voltage-dependent Ca^"*" current /ca,L, and 
subsequently evokes muscle contraction [il,[lB|- A 10% 
stretch is applied to the cell as the stimulation, where the 
strain is defined as the relative change of cell length: 
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FIG. 13. Complete model testing with voltage pulses under 
isotonic condition. A. Ca^^ transient elicited by 1.6 s voltage 
pulse from the holding potential -60 mV to mv. B. The 
cell length change in response to the Ca^^ transient. C, D 
and E. Internal length adjustments of the cell. Isotonic force: 
F =0.4 /iN. 



= ^^£11^ , (30) 
to 

where is the cell length at which the force is zero. 
Figs. [HR-E show the Ca^~''transient and the mechanical 
responses. The underlying electrophysiological events are 
shown in Figs. [TIF-H. An initial increase in force occurs 
immediately after the application of mechanical stretch 
(Fig. [MB). It is due to the quick length (Is) increase of 
the series viscoelastic element (Panel C) and the instan- 
taneous change in length of cross-bridge Ix (Panel E). 
Force drops quickly due to the sliding mechanism of the 
active element (Panel D) which reflects the movement of 
cross-bridges along the actin filament (described as the 
la term in (Eq. This leads to the relaxation of 

the spring element on cross- bridges and series element, 
however, force is still maintained at a lower level until 
Ca^"^ is removed. Membrane depolarization induced by 
stretch-sensitive current Im (Panel G) activates the L- 
type Ca^^ current /ca,L (Panel H), which contributes to 
the intracellular [Ca^^] ;. 



IV. DISCUSSION 

We have taken an integrative modeling approach to the 
characterization of the electrical, chemical, and mechani- 
cal behavior of the single cerebrovascular smooth muscle 
cell. The model takes into account membrane activation 
by either electrical current pulses or mechanical stretch, 
and myoplasmic [Ca^^] regulation, contractile kinetics of 
the actin-myosin interaction, force generation and muscle 
mechanics. 



A. Membrane electrophysiology and fluid 
compartment 

The membrane electrophysiology and fluid compart- 
ment models provide very realistic representations of the 
intracellular Ca^''' transient and hence the input to the 
contractile model. Characteristics of membrane ionic 
currents, especially the voltage-dependent calcium cur- 
rent (/ca,L), were determined using voltage-clamp data 
[O) Uli (see Figs. [3] and |4|). A four-state kinetic model 
was used to mimic Ca^^-induced calcium release (CICR) 
from the sarcoplasmic reticulum (SR), and measured 
voltage-clamp induced [Ca^~'']i transients [13, H^l were 
used to constrain the choice of parameters for the fluid 
compartment and SR models (Fig. [5]). 

B. Myosin phosphorylation and latch bridge 

A combination of hypotheses and measured data were 
used to form a model of active force generation and 
mechanical behavior of the cell. We assumed that the 
Ca'^^ -^-CaCM— s-MLCK pathway is very important in 
myosin phosphorylation and cross-bridge attachment in 
smooth muscle. A modified version of the multi-state 
Hai-Murphy (H-M) model 0, Q was employed to de- 
scribe the kinetics of myosin phosphorylation and cross- 
bridge attachment (Fig. [5]). Specifically, we used the 
activated Ca^~''-calmodulin complex as a trigger for the 
myosin phosphorylation (47j . and made the rate con- 
stants Ki and Kg in H-M model functions of CaCM 
[Eq. (|^^ ]. This implies that in our model, the contrac- 
tile mechanism is tightly coupled with Ca^"*" regulation. 
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FIG. 14. Complete model testing under strain manipulation. 
10% stretch applied to the cell is completed during 100 ms 
period and is held for 1 second before being brought back to 
control level (also in 100 ms). A. Ca^"*" transient in response 
to the strain change. B. Active force (force sustained in active 
branch) response. C, D and E. Internal length adjustment of 
the cell. F. Membrane potential in response to the strain 
change. G and H. Membrane currents elicited by stretch and 
membrane depolarization (/m and Ica.L, respectively). The 
control cell length is 123 ^m. Note that the time scale used 
in Panels F to H is finer than that in Panels A to E. 



Other models, however, consider unphosphorylated cross 
bridge cycling. The cooperative model is an example (48l - 
[50| , whereby it is assumed that nonphosphorylated cross 
bridges are activated, attach to actin and cycle under the 
regulatory influence of a small population of phosphory- 
lated cross bridges. Rembold and Murphy [51j have eval- 
uated multi-state models of latch bridge formation based 
on dephosphorylation and/or cooperativity- regulated at- 
tachment. They consider several four-state cross-bridge 
models that have different kinetic schemes, all of which 
can predict contractile behavior in arterial smooth mus- 
cle. Thus, the process of latch bridge formation is uncer- 
tain and can be either formed by means of dephospho- 
rylation or cooperativity- regulated attachment [5l| . Due 
to lack of quantitative data, we have based our model 
structure on the original kinetic model by Hai and Mur- 
phy B 13 1 and have used the same assumption regarding 
the non-cycling nature of latch-bridges. 



C. Cell mechanics 

We have also developed a model of active force gener- 
ation, and have incorporated it into a mechanical model 
that describes the mechanical coupling of myofilaments 
to so called dense region located among the myofilaments 
and ultimately to the cell wall. Our basic mechanical 
model of the cell is a modified Hill model, which has 
a structure that is similar to that developed by Gestre- 
lius et al. Q. However, it also contains a new descrip- 
tion of the coupling between the contractile kinetics and 
the smooth muscle cell mechanics. Cross-bridge elastic- 
ity and mechanical cycling, as well as, the sliding of cross 
bridges along thin filament are considered, and these pro- 
cesses are modeled as being dependent upon the myosin 
phosphorylation and cross bridge formation kinetics. The 
single cell studies of Harris and Warshaw j^, 01 provide 
length:force (L:F) and force :velovity (F:V) data for vali- 
dation of our mechanical model, and our model provides 
close fits to the measured data (Figs . [TUl and [TT|) . 

D. Model testing and simulations 

Subsystem testing of the integrated model shows: (1) 
the electrochemical model (membrane-fluid compartment 
model) not only provides very good flts voltage-clamp 
data, but to the induced Ca^"*" transients as well; and 
(2) the chemomechnical model (myosin phosphoryla- 
tion/force development/cell mechanics model) provides 
close fits to L:F and F:V curves obtained from single 
smooth muscle cells. The complete integrated cell model 
was tested by predicting the cell response to both volt- 
age pulse stimulation and strain pulses applied to the cell 
membrane with the muscle held under isometric condi- 
tions (Figs. fT2l and fT4|) . Additional tests show that the 
complete cell model can provide simulation of contraction 
of smooth muscle cell under sustained isotonic tension, to 
lift a weight W, through a distance (Fig. [T31). The com- 
ponent isometric and isotonic phases of this response, 
mimic closely what is presently known of smooth muscle 
mechanics. Our model can also predict contractile events 
that are difficult to measure, including: (1) the transition 
occurring between phosphorylated cross-bridges and the 
latch-bridge state during the response; (2) the internal 
length adjustments of mechanical elements within the 
smooth muscle cell, and (3) the integrated responses of 
several different functional components. 

E. Model limitations 

All models have limitations and ours is no exception. 

(1) Our focus in the current study is the development of 
a model for the single vascular smooth muscle cell. This is 
a large task in itself, but the cell model can be integrated 
into larger models of tissue and vessel in which it resides. 
Large quantities of experimental data are available at 
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these more macroscopic levels. The present study should 
be viewed only as a first step in the development of an 
integrative model of the vessel wall, which would neces- 
sarily include functional inter-cellular communication via 
signaling pathways {e.g., nitric oxide, endothelial-smooth 
muscle gap junctions, etc.). Extensions to tissue and 
vessel levels would require additional data regarding the 
connectivity and inter-cellular communication patterns 
involved at these levels. 

(2) Although experimental studies suggest that a slow 
inactivation process is responsible for the reduction of the 
open probability of L-type Ca^"*" channels upon mem- 
brane depolarization 11], quantitative information on 
the voltage dependent time constant and steady-state of 
^Ca.L inactivation is lacking. There is no suitable data 
to quantitatively identify the voltage-dependence of the 
inactivation kinetics. 

(3) The importance of calcium-induced Ca^^ release 
(CICR) in cerebrovascular smooth muscle is in dispute. 
Some groups consider it is important, whereas others dis- 
agree. We have provided a model that exhibits CICR by 
the SR, since CICR may be important in other types 
of smooth muscle cell. The peak magnitude of CICR 
is 10-15% of the peak of Ca^+ transient (Fig. ^. Cer- 
tain experimental studies suggest that: (a) local Ca^"*" 
release through ryanodine-sensitive SR channels (called 
"calcium sparks") are responsible for membrane repo- 



larization via a mechanism whereby the released Ca^"*" 
activates Ca^'''-activated K"*" channels (/K,Ca), which hy- 
perpolarized the cell membrane and limit Ca^^ influx 
In, [13, m, m. These studies also indicate that CICR 
does not contribute to the global [Ca^'^Ji concentration 
in cytosol, and therefore does not directly influence the 
contractile mechanism via the Ca^+-CaCM-MLCK path- 
way; and (b) IPa-mediated Ca^"*" release from SR is al- 
ternative pathway for [Ca^^J^ regulation, which has been 
neglected in the present study, but may need to be con- 
sidered (Hill. 

(4) In our model, we assume that Ca^+-CaCM-MLCK 
activated myosin phosphorylation is the major pathway 
for the contractile kinetics, but alternative pathways may 
exist. Observations indicate that: (i) activation of pro- 
tein kinase C induce slow sustained contractions in vas- 



cular smooth muscle [lj|, and (ii) Rho-kinase sensitizes 
the contractile machinery by inhibiting MLCK phos- 
phatase [5^. Additional data that helps to clarify the 
relative importance of these various pathways, would be 
most important to further clarification of our model de- 
sign. 
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VI. APPENDIX - MATHEMATICAL 
EQUATIONS FOR THE ELECTROCHEMICAL 
MODEL 

Mathematical equations for the model and related pa- 
rameter values are listed in this appendix. 
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TABLE I. Parameters for electro-chemical model 







Numerical value 


"D 

K 


ideal gas constant 


oo4i.U mJ*mol Iv 


1 


absolute temperature 




r 


Faraday's constant 


9D48/.U O-mol 


L'm 


membrane capacitance 


0.03 nF 


5'b,Na 


maximum background Na"*" current conductance 


U.Ui no 


5'b,Ca 


maximum background Ca*^"*" current conductance 


U.Uiz no 


9h,K 


maximum background K"*" current conductance 


U.Ui no 




inward rectifier constant 


U.i4o 


riKi 


inward rectifier constant 


0.5 


pKCa 


maximum Ca^'''-activated K"*" current conductance 


0.5 no 


*^l/2,KCa 


half-activation potential for /K,Ca 


on 

zU.o mv 


r 

JCaP 


ma^ximum Ca'^^pump current 


8.3 pA 


f 

JNaK 


maximum Na~'~-K~'~ current 


/.O pA 


Am,K 


concentration of half-activation for K'^ 


i.U mlVi 




concentration of half-activation for Na"*" 


11. U miVl 


KNaCa 


intrinsic current density of Na^-Ca'^^ exchanger 


U.UUU5 pA miVl 


ONaCa 


scaling factor for Na'''-Ca^''' exchanger 


U.UUUo mJVl 


7 


partition coefficient for Na^-Ca^"*" exchanger 


0.45 


VOlj 


volume of intracellular space 


1.0 pi 


VOlCa 


intracellular volume available to free Ca^"*" 


0. ( pi 


[Oa 


extracellular calcium concentration 


AM mM 




extracellular sodium concentration 


14U.U mlVl 




extracellular potassium concentration 


o.U mlVL 


[ocmJ 


total concentration of calmodulin sites for Ca^^ 


U.l mlVi 


[Bf] 


total concentration of other buffer sites for Ca'^"'" 


0.4 mM 


ki 


Ca^"'"-calmodulin association rate constant 


200 mM-^s"^ 


k-i 


Ca^"'"-calmodulin dissociation rate constant 


0.052 s"^ 


kd 


Ca^^-nonspecific buffer association rate constant 


200 mM^^s"^ 


k-d 


Ca^"*" -nonspecific buffer dissociation rate constant 


0.026 s"^ 



TABLE 11. Cytosolic material balance and calcium buflfering 

dlNa"*"]; _ _ 3/NaK -|- 3/NaCa + .^B.Na + ^M.Na 

di ~ F ■ voli 

d[K+]i _ Ik + JcaK + Im + Ib,k - 2JNaK + Tm,k 
dt_ ~ F ■ voh 

d[Ca^''']i _ Jca.L — 27NaCa + ^CaP + JB,Ca + JlvI.Ca + lup — Irel ^['S'cm] dBp 

dt ~ 2F ■ volca. dt dt 



TABLE III. L-type calcium current 

ICa,L = gCal^dLfhiVm - Eca.,\.) QCah = 1.4151 nS 

ddL dh — dh J 1-0 

dL = 



dt Td^ 1 + e-(^-+l-878)/7.6704 

Td^ = 2.8928e-«^"'+*«3**)/i'-3**-')' + 2.4323 

dff h-ff f 1-0 



■' ]^ _|_ e(Vm+29.3188)/1.5389 

Tf f = 295.5937e~»'^"'~^-^^^^)/"^-^^^>^ 23.1907 /l = 0.74// 0.26 
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TABLE IV. Ca^"*" -activated potassium current 

-^K.Ca = ffKCa-PK,Ca(Kn — Ek) 

_ LO - _ - _ - 

P° " 1 + e-(VV„-Vi/2,KCa)/21™ Pf-Ps- Po 

Vl/2,KCa = -45.01ogio([Ca2+]<) - 198.55 

dpf ^ Pf -Pf dps ^ Ps-Ps 

dt Tpf dt Tn, 

Tpf = 0.5 ms 

PKCa = 0.65p/+0.35p, 



TABLE V. Delayed rectifier 



Ik^ gkPk^iVn^- Ek) 5fe=9.8325 nS 

dpi ^ Pk^Pi ^ 210.9873e-«^™+2" ='=«5)/i'^5-^5''2)' - 20.5866 

dt Tpi 



dP2 ^ Pk -P2 
dt Tno 

1.0 



Pk = fir— T7k ^1/2 = -1-77 mV, A; = 14.52 mV 

Pk = 0.58pi + 0.42p2 



TABLE VI. Inward rectifier 



Vi/2,Ki = 25.19 logio [K+]i - 112.29 



TABLE VII. Pump and exchanger currents 



~ ; CSCKT 

JCaP ~ JCaP 



^^NaK = ^NaK 



CaCM + iCm.CaCM 

[K+]o ^ ( \ fVm + 150 



[K+]o + Km,K ) \ [Na+]^5 + K\^^^^ )\Vm^ 200 
[Na+]?[Ca2+]„$p, - [Na+]2[Ca2+],<I> 



/NaCa - fcNaCa ^ ^ ^^^^^^([Na+T? [Ca^+]» +lNa+]f [Ca^+]o) 



TABLE VIII. Background currents 



lb = /b,Na + /b,K + /b,Ca 

A rv F ^ F _iir [Na+]^ 

[JNa Ji 

/b,K=ffb,K(V^,.-.Bif) SK = ^ln.& 



.^b.Ca = 5b,Ca(Vm — Ecb) Ecb. = 7^ In 2+ 



2F [Ca^+], 
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TABLE IX. Parameters for the multi-state ryanodine receptor (RyR) model 



Parameter 


Description 


Numerical value 




activation rate constant 


2500.0 mM-^-ms"^ 


K,2 


inactivation rate constant 


1.05 mM-i-ms"^ 




unbinding rate constant from activation 


0.0076 ms"^ 


K-r2 


unbinding rate constant from inactivation 


0.084 ms"^ 


I up 


maximum SR uptaking current 


200 pA 


-^m,up 


Michaelis-Menten constant of SR calcium pump 


0.08 mM 


rtr 


time constant of the internal diffusion 


1000.0 ms 


Trel 


time constant of the diffusion from SR release compartment 


0.0333 ms 


vol„ 


volume of uptake compartment 


0.07 pi 


voir 


volume of release compartment 


0.007 pi 


AmPK 


whole-cell K"*" permeability of mechanical-sensitive channel 


1.821e-|-06 L-ms"^ 



TABLE X. Parameters for mechanical model 



Parameter 


Description 


Numerical Value 


hp 


parallel element stiffness constant 


0.1 (uN 




phosphorylated cross-bridge stiffness constant 


12.5 nN/nm 


fex2 


latch bridge stiffness constant 


8.8 ^iN/^im 


kg 


series element stiffness constant 


0.2 /iN 




viscosity of series viscocritical element 


0.01 /iN-ms//xm 


^0 


cell length at zero passive force 


40 /im 


ho 


length of series viscoelastic element at zero force 


30 /Ltm 


lopt 


optimal length of active contractile element 


100 urn 


/amp 


friction constant for phosphorylated cross-bridges 


1.3 /iN-ms//im 


/am 


friction constant for latch bridges 


85.5 /nN-ms/yum 


Vx 


cross-bridge cycling velocity 


5.0 /im/ms 




length modulation constant 


7.5 




constant for parallel clement 


0.1 


as 


constant for series elastic element 


4.5 



TABLE XL Parameters for myosin kinetics model 



Parameter 


Description 


Numerical value 


K2 


myosin dephosphorylation rate constant 


0.4 s"^ 


Ks 


cross-bridge formation rate constant 


1.8 s"^ 


Ki 


cross-bridge detachment rate constant 


0.1 s-i 


Ks 


myosin dephosphorylation rate constant 


0.4 s"^ 




latch state to free myosin detachment rate constant 


0.045 s"^ 


J^CaCM 


half activation constant for myosin phosphorylation 


1.78e-7 
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TABLE XII. Initial conditions for state variables 



-56.1275 mV 

li=100.92 nM 



Vm 

[Ca 

[Na+]i = 5.06 mM 
[K+]i = 91.03 mM 
dL = 0.00046 
// = 1.000 
fs = 1.000 
pi = 0.0159 



P2=0.0159 

p/=0.00032 
Ps=0.00032 
[S'cm]=0.072 mM 
[Bf1=0.336 iriM 
M = 0.6099 
Mp = 0.0476 
AMp = 0.0627 



7?oi =0.9955 
7?io=0.0033 
i?ii=4.0e-06 
[Ca^+]up = 0.5687 mM 



,2+1 



is=37.05 /um 
=89.60 Aim 



=0.5503 mM 



